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Abstract. Dynamical systems studies of differential equations often focus on the behavior of so- 
lutions near critical points and on invariant manifolds, to elucidate the organization of the associated 
flow. In addition, effective methods, such as the use of Poincare maps and phase resetting curves, 
have been developed for the study of periodic orbits. However, the analysis of transient dynamics 
associated with solutions on their way to an attracting fixed point has not received much rigorous 
attention. This paper introduces methods for the study of such transient dynamics. In particular, we 
focus on the analysis of whether one component of a solution to a system of differential equations can 
overtake the corresponding component of a reference solution, given that both solutions approach 
the same stable node. We call this phenomenon tolerance, which derives from a certain biological 
effect. Here, we establish certain general conditions, based on the initial conditions associated with 
the two solutions and the properties of the vector field, that guarantee that tolerance does or does 
not occur in two-dimensional systems. We illustrate these conditions in particular examples, and we 
derive and demonstrate additional techniques that can be used on a case by case basis to check for 
tolerance. Finally, we give a full rigorous analysis of tolerance in two-dimensional linear systems. 
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1. Introduction. Relative to asymptotic behavior, transients have received ht- 
tle attention in the study of nonlinear dynamical systems. For example, how the rate 
of approach to a stable fixed point, away from the asymptotic limit, is affected by 
the choice of initial conditions within the basin of attraction of that fixed point has 
not to our knowledge been well characterized. In this work, we consider a comparison 
of the transient dynamics of pairs of trajectories with similar asymptotic behaviors. 
The motivation for this work arises from a biological phenomenon known as tolerance, 
which refers to a reduction in the effect induced by the application of a substance, 
due to an earlier exposure to that substance. For example, administration of a toxin 
to rodents, at a given reference dose, induces a reproducible acute infiammatory re- 
sponse featuring a rise in a variety of immune system elements followed by a return 
to near-baseline conditions [U HI [TTJ [T3] • If a small pre-conditioning dose of the toxin 
is given to an animal prior to the reference dose then the activation of immune agents 
by the reference dose is attenuated. This phenomenon is called tolerance. 

A previous study [5] analyzed tolerance in the context of a four dimensional ordi- 
nary differential equation (ODE) model of the acute infiammatory response. Within 
the four dimensional ODE model, the origin represents a healthy equilibrium state, 
and the abrupt administration of a toxin is represented by a jump of a trajectory to 
another point in phase space. Thus, starting from a given initial condition, tolerance 
occurs precisely when the sequence of a pre-conditioning dose, a period of ensuing 
fiow, and a subsequent reference dose leads to a trajectory position that is different 
from the one attained by direct administration of the reference dose, and from which 
a lower level of activated immune agents ensues. From the observation of tolerance 
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in the acute inflammatory response model, we reasoned that similar tolerance effects 
should be a general feature of trajectories generated from different initial conditions 
by a dynamical system with negative feedback. Little analysis has been done on tran- 
sient effects such as tolerance, compared to the major emphasis in dynamical systems 
research on invariant manifolds and other structures derived from asymptotic and 
local calculations f8^, TT . 

Our goal in this work is to provide a framework for the study of tolerance in 
ODE systems. Specifically, we focus on trajectories converging to an asymptotically 
stable node. Overall, we are interested in necessary and sufficient conditions for 
tolerance, as we formally define it in Section 2. In a one-dimensional or scalar ODE, 
uniqueness of solutions prevents tolerance from occurring. Thus, we examine tolerance 
in two-dimensional ODE systems, using geometrical approaches. The general two- 
dimensional nonlinear case, which is treated in Section 3, poses challenges, since 
exact analytical solutions are generally not available. However, through the use of 
isoclines and the concept of inhibition, we give some general results on conditions 
when tolerance can or cannot occur and we develop an approach to the derivation 
of more precise results for particular models. Specific examples axe used here to 
illustrate this approach. In Section 4, we take advantage of analytical solutions to 
provide a complete analysis of tolerance in two-dimensional linear systems. We finish 
with conclusions and a brief discussion of related work in Section 5. 

2. Preliminaries. 

2.1. Definitions and assumptions. In this section we present our assumptions 
and give the precise mathematical definition that we use for tolerance. Consider the 
autonomous ODE system 

i = f{x,y) (2.1) 
y = 9{x,y), 

where € M, and f,g are locally Lipschitz. 



(Al) Assume that there exists a stable fixed point of (2.1 ), the eigenvalues of which 
are real and negative (to eliminate spirals and centers). Without loss of 
generality, we will take (0, 0) as the given stable fixed point of ( |2.1[ ). 

Let be the basin of attraction of (0, 0) in the first quadrant, l^f [0, oo) x 

[0,oo): 

rfo.o) = n {(x, y)\{x, y)-t^{0,0) ast^ c»}, 



where the notation (x, y) ■ t is the image of the point (x, y) under the flow of (2.1 1 for 
time t. The set of points, {{x,y) ■ t\t > 0}, is the solution curve or trajectory of the 
initial value problem with initial value {x,y). This set is also referred to as the graph 
of the solution. 

Let 4){t) = (0i(i), (/)2(^)) and -ipit) = (V'i(t), V'2(i)) be two solutions to the initial 



value problem of (2.1) with initial values 

(/.(O) = (x„y,), X, > 0, > (2.2) 

and 

■0(0) = (xp, yp),Xp >0,yp> 0. (2.3) 
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(A2) Assume that both components of ip{t) and ip{t) are nonnegative for all i > 
and that (x^jUr) and {xp,yp) G T^p^. 

(A3) Assume that Xr and Xp are chosen such that Xp > Xr- 

Definition 2.1. Define (j)(t) as the reference (R) trajectory or solution. 

Definition 2.2. Define ip{t) as the pre-conditioned or perturbed (P) trajectory 
or solution. 

Essentially, we are interested in determining whether or not there exists a time 
when the first component of a P trajectory overtakes that of an R trajectory, given 
that it was initially behind, as they approach the origin. Our ensuing discussion would 
apply equally if we considered the second component instead of the first. 



Definition 2.3. The system (2.1) is said to exhibit tolerance for {{xr, Vr), {xp, yp)) 
if there exists r > such that ?/'i(r) < 0i(r). 



Definition 2.4. If ipiit) > (l)i{t) for all t e [0,oo), then (2.1) does not exhibit 
tolerance for {{xr,yr), {xp,yp)) . 



Remark 1. We will also use the terminology that {xp,yp) or ip produces (or 
does not produce) tolerance in (2.1 ) with respect to {x.^ yr) or (p to mean that Defini- 
tion 2.3 (Definition 2.4) holds. Figure 2.1 illustrates definitions \2.3\ and 2.4 with time 
courses of the first component of solutions (f>{t) and ^p{t) for a given {{xr, yr), {xp, yp)). 





Fig. 2.1. Illustration of Definitions \273\ and \2.4\ Left (Right) panel: Time course of the first 
component of a pre-conditioned (P) solution, 1/1(4), with initial condition (xp,yp), which produces 
(does not produce) tolerance with respect to the reference (R) solution, <j>(t), with initial condition 



Remark 2. Under (A3), V'(O) {xp,yp) e [^^,00) x [0,cx)); that is, the initial 
value for the P solution could lie at any point on or to the right of the line x = x^ in 
the first quadrant. Correspondingly, we define ^^qq) to he the basin of attraction of 
(0,0) in [xr,oo) X [0,00) C M^+z 

^Fo.o) = ^(0,0) [2^r,oo) X [0,00). 
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Remark 3. The above definitions of tolerance are related to the biological setting 
that motivated this study through the interpretation of the P trajectory. Consider a 



non-negative pre-conditioning solution pit) — (pi(t), p2(t)) of (2.1) with initial value 



p(0) = {Xp,yp), < Xp < Xr, < Vp < Vr- 

We then interpret the pre-conditioned solution — {tpi(t) , tp2(t)) as the solution of 
(2.1) with initial value 



V'(O) = {xp{s),yp{s)) = p{s) + ixh,yh) for some < s < oo, 



(2.4) 



where {xh,yh) G R^"*". If{xh,yh) = {xr,yr), which is typical for inflammation exper- 
iments, then for fixed (t>{0) = {xr,yr) and p(0) = {xp,yp), every s defines a unique 
initial value for ijj that satisfies (A3), namely {xp{s),yp{s)) as defined in equation 
[2.4-). Thus, for a continuum of s values ranging from to oo, a curve of possible 
{xp, yp) values is formed, and it is of biological interest to know which of these {xp, yp) 
lead to tolerance. 



2.2. Properties of tolerance. Definition |2.3| refei s only to the presence of toler- 
ance at one time point t > such that iPi{t) < 0i(t). However, continuity arguments 
can extend this window from a single time point to an open interval, (ti,t2), around 



T, with V''i(^i) = 01 (^i)- This observation is stated formally in Proposition 2.5 below 



and will be important in Section [3] Figure 2.2 illustrates Proposition 2.5 with time 
courses of relevant solutions. 



Proposition 2.5. Assume {Al), {A2), and {A3). If (2.1\ exhibits tolerance for 
{{xr,yr), ixp,yp)) at T > 0, then there exists an open neighborhood (^1,^2) around 
T such that "01 (i) < 4>i{t) for every t € (^1,^2) o,nd = (/)i(ti). Furthermore, 

fmh))<f{<t>{t,)). 




Fig. 2.2. Time courses illustrating Proposition \2. 5\ Note that in this example, t2 could be 
chosen to be any t > t. 

The window of tolerance can also be extended with respect to (xr, yr) and (xp, yp). 



Proposition 2.6. Assume (Al), {A2), and {A3). If (2.1\ exhibits tolerance for 



{{xr,yr), {xp,yp)), then there exists an open ball, Br, of radius r around (xr,yr) such 
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that if {xk,yk) G Br{{xr,yr)) ^'^^Q^y then there exists a corresponding time tk > 
such that tolerance is exhibited for {{xk,yk), {^ptVp)) ■ 



Proposition 2.7. Assume (Al), {A2), and (AS). If (2.1) exhibits tolerance 
for given {{xr,yr)i{xp,yp)) , then there exists an open ball, Bf, of radius r around 
{xp,yp) such that if {xk,yk) G Bf{{xp,yp))C\T^^^^y then there exists a corresponding 

time tk > such that tolerance is exhibited for {(xr,yr), {xk,yk))- 



Propositions 2.6 and 2.7 are easily proved by noting that solutions of (2.11 are 
continuous and depend continuously on initial conditions. Each time tk or tk can also 
be extended to an interval of times for which tolerance occurs, by Proposition |2.5[ 



3. Conditions for the existence of tolerance. In this section, we progres- 
sively build up a collection of ideas that are useful for determining the set of initial 
conditions for P for which tolerance can be guaranteed to occur or not to occur. In 
particular, in subsection |3.1[ we present a basic result on a general situation in which 

we introduce some concepts 



tolerance can be guaranteed to occur. In subsection |3. 2 



that are useful for refining the results from subsection 3.1 and we discuss their im- 
mediate consequences for tolerance. We harness these ideas in subsection 3.3 where 
we set up a general approach that can be used to move beyond the results from sub- 
sections |3.1| and |3.2| in particular systems, and we illustrate this approach in several 
examples in subsection |3.4[ 



3.1. Basic conditions. In this subsection, we consider specific conditions on 
the initial values of P and R for which tolerance can or cannot occur. We first con- 
sider conditions in which tolerance can occur when solutions (j)(t) and ip(t) of system 



(2.11, as defined in Section 2.1 are subsets of the same solution curve. 



Proposition 3.1. Assume (Al), {A2), and {A3). Given {{xr,yr), {^piVp)) ! o*" 
sume (j)i(t) and ip ijt) monotonically as t —> oo. // there exists t > such that 



4i{~i) — (xp,yp), (2.1) does not exhibit tolerance for {(^Xrjyr), (^ptVp)) ■ 



This proposition follows immediately from the group property of flows and is the 
reason why tolerance is ruled out in one dimensional systems. Next, we focus on a 
situation where the reference trajectory is what we call an excitable trajectory as 
represented, for example, in the left panel of Figure |3.1[ We make this precise in 
terms of the graph of given by 

graph{4>) = {{x, y) = {xr,yr) • t : t > 0} , (3.1) 

with the following definition. 

Definition 3.2. Assume that [Al), {A2), and (AS) hold. Fix a positive integer 
n. The trajectory (l){t) is n-excitable if there exist times t^g — 0, , . . . , tg^^.j > 
such that 

(a) (}ii(iei) > Xr for all i > 0, 

(6) g{Mt),Mt)) >0 forte [0,te.„_J, and 

r f{Mt),Mt)) >0, ie [te„,te,) and {te,,,t,,^^,), i G {1, 2, . . . , n - 1}, 

\ f{Mt),Mt)) < 0, t e (te2,+i,ie,(,+i)), i e {0,l,...,n-2}, ort> te^^_^. 
The trajectory (j)(t) is excitable if it is l-excitable. 
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Excitable trajectories are common in various biological models. In the context 
of acute inflammation, an excitable trajectory represents the initial activation of the 
immune system by a stimulus followed by a relaxation to a stable baseline state. 



Remark 4. Condition (b) on g in Definition 3.2 is not necessary for our ap- 
proach, hut this assumption clarifies the presentation to follow. 

Below, we define a set T such that tolerance with respect to (x^, Ur) occurs when- 
ever {xp, Up) € T, when (f)(t) is an n-excitable trajectory. 

Definition 3.3. For an n-excitable trajectory cj), define tr > to be the first 
positive time where 0i(tr) = Xr, which exists since <f> is n-excitable and continuous 
{and by {Al) and (^2)). Note also that (pi(t) > 4>i(tr) = Xr for all t € (0,fr) by 
definition of an n-excitable trajectory. 

Definition 3.4. Now, in terms oftr, define G to be the set of points {x,y) ^ 
{xr, Vr) on the graph of (f> for t € (0, tr] : 

G = {{x, y)\{x, y) = c^{t) for t e (0, tr]} . (3.2) 



Definition 3.5. Assume that <f) is an n-excitable trajectory. Define L to be the 
line segment L — {x : x — Xr,y € {yr , (l)2{tr)]} and define the region S {see Figure 



3.1] as the union of L and the interior of the region bounded by G and L. 

Definition 3.6. Define T as the union of G and S as defined above, 

T=GUS. (3.3) 



Definition 3.7. Define M — maxf>o{0i(t)}, which exists by {Al), {A2), and 
the continuity of (f). Let t^ > {tn > 0) be the minimal (maximal) positive time such 
that (j)i{t) = M. 

Proposition 3.8. Let (j){Q) — {xr,yr) and let {xp,yp) be given. Suppose that 
{Al), {A2), and {A3) hold and that cj) is an n-excitable trajectory. Under these condi- 
tions, T is a non-empty set. Moreover, if{xp,yp) S T, then ^2.1\ will exhibit tolerance 
for {{xr,yr), {xp,yp)). 



Proof. By the assumptions, a region T — GUS as defined above exists. We divide 
the proof into two parts since T is defined as the union of two sets. 

Part 1: Suppose ip{0) = {xp,yp) G G. This implies that ip{Q) = {xp,yp) = 4>{t), 
for some r > 0. Again, <^i(t) < M for all nonnegative t > Im. It follows that V'i(^m) = 
01 (^M + t) < M — 4>i{tM). Thus, (2.1 ) exhibits tolerance for {{xr,yr)i {xp, yp)) G G 
at time ^m- 

Part 2: Suppose (xp, yp) G S. We first consider the case where Xp > Xr and define 
tp = minf>o{i : i'i{t) = Xr} such that ')p{t) G S for all t G [0,tp]. If tp > tr then since 
tr > tM > tm, tm G {0,tp). Hcncc, ijji{t,n) < M = 01 (im) and tolerance is exhibited 
at tm. Now, if < < tr, then it is possible that 'ipi{tm) > M (see bottom panel 
of Figure 3.11. However, from the definition of i,., (f>i{tp) > 4'i{tr) = Xr = V'i(^p) 
and tolerance is exhibited at tp. Now, consider the special case that Xp — Xr. If 
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f{xp, Hp) > then one of the above two cases holds. If f{xp, yp) < 0, then there exists 
e > such that ipi{e) < Xj. and 0i(e) > x^.- Thus, <j){e) > '0(e) and tolerance occurs 
at e. □ 



Figure 3.1 illustrates Proposition 3.8 in both phase space (left panel) and with 



time courses (right panel). Notice that if we consider the special case when (xr,yr) 
of an n-excitable trajectory is on the x-axis, then uniqueness of solutions is sufficient 
to guarantee tolerance. 



/ \ 












(X 







Fig. 3.1. Illustration of Proposition [3. 8\ in the case that (p is n-excitable. P trajectories with 
initial conditions in region S exhibit tolerance. Left Panel: A 2-excitable R trajectory, initial 
condition, {xr,yr) (black) and two example P trajectories, ip{t), initial condition {xp,yp) £ S (red). 
The maximum value in the x-direction for <j>{t) is marked with a vertical blue line and denoted by 
M . Right Panel: Time courses of both <j>\{t) (black) and (red). Time tp is where jpi first takes 

on the value Xr and t^ is time when <pi (t) last attains its maximal value. 

If more constraints are imposed on the vector field / then the region that guar- 
antees tolerance can be immediately expanded to include the strip above T in T^^ . 
To be precise, we introduce the following definition. 

Definition 3.9. Define f by the set 

f = {{xr, M) X (02(iAf), oo) \ T) n Vl^^^y (3.4) 

Proposition 3.10. Assume {Al), {A2), {A3), and that 4> is a n n-e xcitable tra- 



jectory with 0(0) = {xr,yr)- If f inT, then for {xp,yp) € T, (2.1\ will exhibit 
tolerance. 
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Proof. For {xp,yp) G T and / < 0, it follows from the assumptions that %IJi{t) < 
< 0i(^A/) for t > 0. Thus, 01 (iA/) > '0i(^a/)- Hence, (2.11 exhibits tolerance for 
{{xr,yr), {xp,Vp) e T > at time tM- □ 



3.2. Isoclines and Inhibition. In the previous section we found generic con- 
ditions under which tolerance would occur. However, the initial conditions resulting 
in tolerance were confined to a small region of the available basin of attraction. Nu- 
merical experiments in various examples suggest that the region for tolerance is often 
larger. Here, we introduce new concepts that enable us to expand the regions on 
which we can show that tolerance is possible or guaranteed. 



Consider the ODE (2.1) and assume (Al), (A2), and (A3) hold. 



Definition 3.11. The a; -isoclines o/ (2.1\ are the family of curves {or level 



sets), parametrized by a parameter C G M, each defined by f{x,y) = C. 



A nuUcline, for instance, is an isocline for which C = 0. The vector field points 
in the positive (negative) s-direction when C is positive (negative). 

Remark 5. We may define y-isoclines analogously to x-isoclines. Since we do not 
consider these, we will drop the x- and just use isocline to refer to the x-isoclines here. 

We now introduce the concept of inhibition. Inhibition is a widely used term, 
especially in the context of mathematical models of biological systems, for the sup- 
pression of one quantity by another. However, the use of this term, while intuitive 
and heuristically understood, is not always mathematically precise. Hence, we give a 
precise definition of inhibition. Subsequently, we prove two results relating to inhibi- 
tion and tolerance. 



Definition 3.12. Given ft C R-^+ , y inhibits x in ft, and ft is a region of inhi- 



bition/or (2.1), if f{x,y) is a monotone decreasing function of y in f2. 



Remark 6. Note that the sign of f{x,y) is not specified in Definition 3.12 
when y inhibits x, it may either slow the growth of x or speed up its decay. 



Thus, 



A key first observation that follows from the definition of inhibition is that there 
is always the possibility of tolerance when y inhibits x, as long as the perturbed tra- 
jectory samples larger y values than the reference trajectory. We now formalize this 
observation by stating two further definitions and proving two preliminary results, 
which establish the necessity of a region of inhibition and of certain relative positions 
of the perturbed and reference trajectories, respectively, for tolerance to exist. 

Definition 3.13. The graph oftp is bounded below by the graph of4> if 4'2{si) < 
ip2is2) whenever 0i(si) = ^/'i(s2) for any si, S2 > 0, not necessarily equal. For brevity, 
we say is bounded below by 0. 



Proposition 3.14. Assume that (Al), {A2), and {A3) hold and that tp is 
bounded below by 4>. If (2.1\ exhibits tolerance for a given pair {{XrTyr),{Xp,yp)) , 
then there exist a region of inhibition il and si,S2 G such that ^i{si) = 0i(s2) 
with 7/'(si), 0(52) G 
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Proof. Assume that tolerance exists for {{xr,yr), {xp,yp)) but y does not inhibit 
X in any region fl that contains points (V'i(si), -02(51)) and (0i(s2), 02(52)) where 
V'i(si) — 0i(s2) and si,S2 G M+. Given tolerance, it follows from Proposition 2.5 



that there exists t* such that ipi{t*) — 0i(t*) and ipii^t) < 0i(t) for all t e {t*,t* + S) 
for some S > 0. Thus, f{ijj{t*)) < f{(j){t*)). Since the graph of is bounded be- 
low by the graph of 0, we have that at t* , 'ip2{t*) > 02 (i*)- Our assumption that 
y does not inhibit x in any region il containing the points '0(^*) and 0(i*) implies 
f{ip{t*)) > f{(j){t*)), which is a contradiction. Hence, if ip is bounded below by 0, and 
(2.1 1 exhibits tolerance for {{xr,yr)i{xp,yp)) , there must exist a region of inhibition 
and si,S2 G K^, such that ?Ai(si) = 0i(s2) and V'(si),0(s2) S fi. □ 



Remark 7. A^ote </ia< Propositions \3.r^ and \3.1J\ together imply that for an n- 
excitahle trajectory to exist there must exist a region of inhibition. 

Proposition |3. 14| states that a region of inhibition is necessary for tolerance to oc- 
cur when the P trajectory, V'(^); is bounded below by the R trajectory, 0(i). However, 
for ^ bounded above by 0, inhibition can be a detriment to the presence of tolerance 
under certain conditions. First, we define what it means for ip to be bounded above 
by 0. 



Definition 3.15. The graph ofijj is bounded above by the graph of(f> z/02(si) > 
^^2(32) whenever 0i(si) — "01(52) for any si, S2 > 0, not necessarily equal. For brevity, 
we say ip is bounded above by 0. 



Proposition 3.16. Assume that [Al), {A2), and {A3) hold. For {{xr,yr),{xp,yp)) 
such that Xp > M, if the graph of ip is bounded above by the graph of 4>, and y inhibits 
X in a region f2 such that (j){t), ip{t) C fl for all t > 0, then ^2.1\ cannot exhibit 
tolerance for {{xr,yr),{xp,yp)) ■ 



Proof. The proof is analogous to that of Proposition |3.14 If Xp > maxt>o 0(t) 



then tolerance requires f{'p{t*)) < f{(t>{t*)) for some t* such that 0i(t*) = 0i(t*), but 
this cannot occur in a region where y inhibits x, given that ip is bounded above by 0. □ 



Thus, Proposition |3.16| states that in order for tolerance to be a possibility for 
a P trajectory tp that is bounded above by the R trajectory 0, for initial condition 
V'i(O) > M, it is necessary that there exists at least one pair, si,S2 € such that 
V'i(si) = 01(52) and -0(51), 0(52) do not belong to a region of inhibition. 

Propositions |3.8| |3.10| |3.14[ and |3.16| suggest a strategy for evaluating whether 
or not tolerance may occur in a particular system for given R and P trajectories with 
initial values 0(0) = {xr,yr) and "0(0) = {xp,yp), under assumptions (Al), (A2), and 
(A3). First, if is an n-excitable trajectory, then b y Pr oposition [3^ tolerance occurs 
for all (xp, j/p) e T (see Definition 
then by Proposition 



3.10 



3.6 



and Figure 3.1 ). If in addition / < in T, 



tolerance occurs for all {xp,yp) £ T (see Definition 3.9 1. 
Next, we identify the regions of inhibition for system (2.1 1. If it can be established 
that the trajectory ip emanating from an initial condition (xp,yp) is bounded below 
by but does not pass through a region of inhibition, then tolerance cannot occur 
(see Proposition 3.141. Similarly, if Xp > M, ip is bounded above by 0, and -0,0 are 
contained in a region of inhibition, then tolerance cannot occur (see Proposition 3.16 1. 
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If /j, < on all of M^+, then the possibility of tolerance exists for all {xp,yp) such 
that ip is bounded below by (j). 



3.3. Time interval estimates. To obtain more precise conditions for the exis- 



tence of tolerance, direct estimates regarding specific trajectories of ( 2. 1 ) are necessary. 
Here, we show how to derive estimates for upper and lower bounds on the amount of 
time it takes for the relevant trajectories to reach a specified a;- value Xf that is crossed 
by both trajectories, (j)(t) and ip{t). If an {xp,yp) can be found such that Tp{t) takes 
a shorter time interval to reach xj than (/'(i), then tolerance exists for that {xp,yp). 

Assume that there is a positive integer n for which the graph of (fi can be de- 
composed into a union of n graph segments such that the y component of the graph 
is single valued with respect to x on each. This assumption holds, for example, 
when (p is m-excitable for some m. Let Xi, i G {1, . . . ,n -I- 1} be the n -I- 1 termi- 
nal points of the n segments, defined by xi = x^, Xi — for i = 2, ...,n, 
where = infj^^i-ijt : 02(i)) = 0} with ij, 0, and Xn+i — xj. Let 

t^"*"^ = inft>t^{t : = Xf}. The total time to traverse the trajectory from Xr to 

Xf \s then given by = X]"=i ^^0; where A<*^ = - f^. 

On each graph segment we can express the graph of as a function y — Vi{x), 
where Vi is defined on the interval Xi < x < x^+i, i S {1, . . . , n}. We can compute Ati 
for each segment directly by integrating the first equation of (2.1) along the graph 
segment defined hy y — Vi{x), i.e. x = f{x, Vi{x)), to obtain 



f{U,Vi{u)) 



A similar construction can give t^, with initial x-coordinate Xp. Tolerance then implies 
tijj < tcf,. In general, it is not possible to obtain Vi in closed form, but depending on 
the structure of /, estimates can be made to obtain various bounds for and t^. 

For example, with respect to ( |2.1[ ), consider the family of x-isoclines f{x,y) — C, 
where C G M. Let = sup^^j^i t>+i){|/('/'i(i)j </'2(i))|}, i-e. the largest magnitude iso- 
cline through which the trajectory (j) passes on the segment [x\, x]^^]. Then from (3.5 1 



we obtain > J2'i=i W^^-x\\lc\- Likewise, let c\ = inftg[jj^_jj+i){|/(-0i(t), ■02(i))|}, 
i.e. the smallest magnitude isocline through which the trajectory '0(i) passes on the 
segment {x\,x'^\ yielding < XI^i Ix'^^ - Thus, if 



then < tth, which implies tolerance. 



1=1 



We can use condition (3.6 1 to show, for example, that if ipit) is bounded below 
by an m-excitable trajectory (f>{t), and cjj^t) and '>p{t) both lie in a region of inhibition, 
then the region on which tolerance is guaranteed to occur can be expanded from 



that defined in Proposition 3.10 As an example, suppose that (j){t) is an excitable 
trajectory. We can then divide (f) into two segments. In the first segment (pi{t) and 
4>2(t) are increasing, and in the second (l)i{t) is decreasing. By continuity and (Al), 
4>2(t) must first increase and then decrease on the second segment. The end point 
of the first segment is xm — niaxt>o (/'i(t). Define Xf as the x-value where (t>2it) is 
maximal and let 02 (^) — Vf Sit this point. Since (f>(t) belongs to a region of inhibition. 
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the largest magnitude isocline through which the first segment of (p{t) passes is given 
by cj, = f{xr, Vr) = Cr- On the second segment, the largest magnitude isocline passes 
through (t){t) when (t)2{t) is maximal. Thus c? = |/(a;/,j//)| = Cf > 0. 




iO.O) 



Fig. 3.2. Illustration for time interval estimates. 



Now, using Figure 3.2 as a reference, consider a trajectory ip{t) such that / < 



along the trajectory, so there is only one segment and it is bounded below by the line 
y = ijf. Thus, c\ — > Cf, and tolerance is observed if 



\xf - Xp\ ^\xM - Xr\ , |a:;/-a:;A/| 



Cr 



Cf 



(3.7) 



If we consider an excitable trajectory, then Xp > Xf, xm > Xr, and xm > Xf. Taking 
these inequalities in (3.7 1 gives the tolerance condition 



Xp < XM + 



Ctb — C 



Cf 



(xm - Xf) + -^{XM - Xr) = XM 

O7- 



(3. 



Since > Cf, (|3.8[ ) implies that X]\i > xm , which expands the region obtained from 
Proposition 3.10 We note that C^ is a function of yp, so (3.8 1 defines a region R such 
that if {xp, yp) € R, then tolerance occurs in (2.1 1. 



3.4. Examples. In the examples below, we illustrate the ideas introduced in the 
previous subsection. 

Example 3.17. Consider the system given by 



X 

y 



.fix,y) 
9{x,y) 



l+y 

-T.2 _ S 
X 2 



X 



(3.9) 



Note that (0, 0) is a stable node for (3.9 1. The isoclines for this system are the family 
of curves given by the equation 



-x-C 
x + C 



(3.10) 
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for C e M. Figure 3.3 shows a subset of the isodines for C E [—4.0,50] shown in 
increments of 0.5 for those above the C = isochne and in increments of 1.0 for those 
below. 



30 



Entire first quadrant 
is a region of inhibition 




C<0 



C>0 



1 2 3 4 5 6 7 8 

Fig. 3.3. Isodines for Example \3.11\ defined by Equations \3.1C\ for C S [—4.0, 50]. 



For each C < 0, the corresponding isochne has a local minimum at a; = — 2C and 
a vertical asymptote at x — —C. Direct differentiation of / in (3.9 1 yields fy < 0, or 
equivalently, from (3.10), dy/dC < 0, for all {x,y) in the first quadrant. Thus, the 



entire first quadrant is a region of inhibition. We will consider several different initial 
conditions {xr,yr) for (pit) in this example: 



(a) = (4.0,0.0), 

(b) = (4.0,3.0), 

(c) (x,,y,) = (4.0,10.0). 



For initial condition (a), Figure 3.4 'a) displays the following features: 

{Xr.yr) 



4>{t) is the curve shown in black for initial condition (j){0) = {xr, yr) = (4.0, 0). 
Ri is the green region vmion the boundary of 4>{t) and is defined in the same 
manner as the region T, in Definition |3.6[ Ri is bounded to the left by 
{x = Xr}, in accordance with (A3), and to the right by {x — xm}- 
i?2 (the light blue region) is the strip in r^jJ^Q^ lying above sharing its 
bounds on x. 

R3 (the yellow region) is the complement of Ri U i?2 with respect to F^^ , 
namely R3 1^^ r^o^o)\(i?i U R2). 



Case 1(a): {xp,yp) g By Proposition 3.8 any (xp,yp) € Ri will produce 
tolerance. Furthermore, define G — {(x,y)\(x,y) — graph((j)) n [x,-,oo) x (0, 00)}. By 



Proposition 2.7 for each {xp, yp) £ G, there exists an open ball, Bf., of radius f around 
{xp,yp) such that {xk,yk) € Bf C^^^qq-^ produces tolerance with respect to (4.0,0). 
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(a) V 

Entire first quadrant 










\ C = -4.0 












^^^^^^^ 


Ri 

J 





%,^vH4,0) 



(b) 



Entire first quadrant 
is a region of inhibition 



C<0 
OO 

X 




(c) 



Entire first quadrant 
is a region of inhibition 




Fig. 3.4. Isoclines and various initial values {xr,yr) for Example\3. n\ a, b, and c. 



Case 2(a): {xp,yp) £ i?2- Region i?2 is a region of inhibition in which / < 0. 
Thus, by Proposition 3.10 any {xp,yp) G R2 will produce tolerance. 



Case 3(a): {xp,yp) e R3. In this case, for ipi{t) > AI ^ maxt>o{0i(i)}, ij{t) 
is bounded below by 0(t) and the presence of a region of inhibition makes tolerance 



possible (Proposition 2.7l. For ijji^t) < M, which is possible for small j/p, 'ip{t) will 
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1 2 3 4 5 6 7 8 











Fig. 3.5. Animations for Example \3 . n\ ^a) , showing the presence or absence of tolerance with 
respect to 0(0) = (4, 0) for differing choices ofip{0). 0(0) is denoted by the large red dot and V'(O) is 
denoted by the smaller blue dot. Given (f>{0) = (4.0,0), ^l>{0) = (4.5,5) £ Ri produces tolerance {Top 
Left), "0(0) = (4.5,20) S i?2 produces tolerance (Top Right), i/'(0) = (6,10) £ R3 produces tolerance 
{Bottom Left), and V'(O) = (7, 1) S -R3 does not produce tolerance {Bottom Right). 



eventually be bounded below by (pit) and hence tolerance is again possible. 



Figure |3.5| contains links to four separate animations that illustrate the presence 



or absence of tolerance in Example 3.17 'a) using various choices of ip(0) from the 
different regions shown in Figure 3.4 ^a). Each animation displays both phase space 
trajectories of and ip and time courses of and tpi{t) in a side-by-side compar- 
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If ijr is increased with Xr fixed, the regions i?i and R2 shrink. Finally, when yr 
reaches 3.0, corresponding to initial condition (b), these regions disappear. Figure 



3.4 'b) displays the following features: 

• (f>{t) is the curve shown in black for initial condition 0(0) = {xr,yr) = 
(4.0,3.0). 

• The orange curve, denoted as 0, is the curve of points obtained by integrating 
(j){t) backwards in time from t — to t ^ —1.0, at which time it intersects 
the X-axis at a; ~ 3.4. 

• i? is the yellow region defined to be F^gg^ \ (4.0,3.0). 

For this example, if Xp = = 4.0, then for all {xp,yp) G R, the corresponding graph 
of is or will eventually be bounded below by the graph of 0. Since the graph of "0 lies 



in and M^"*" is a region of inhibition. Proposition 3.14 implies that it is possible 
that tolerance can be exhibited by any [xp, yp) G R, although, as in the previous case, 
tolerance is not guaranteed (see Figure [Oj^b)). 



For yr > 3.0, the situation is qualitatively similar to that shown in Figure 3.4 c) 



for initial condition (c), {xr,yr) = (4.0,10.0). Figure 3.4 'c) displays the following 
features: 

• 4)(t) is the curve shown in black for initial condition 0(0) = {xr,yr) = 
(4.0,10.0). 

• The orange curve, denoted as 0, is the curve of points obtained by integrating 
(p{t) backwards in time from t — to t ^ —.58, at which time it intersects 
the X-axis at a; « 4.0. 

• Ri is the orange region union its boundaries: (1) 0(t) and (2) the line segment 
{{x,y)\x = xr,y e [0,10]}. 

• i?2 is the yellow region defined to be the complement of i?i with respect to 
r^ov namely R2 = Fj"'' . \ 



Case 1(c): {xp,yp) £ Ri. Using Proposition 3.16 and Proposition 3.1 {xp,yp) G 



i?i cannot produce tolerance with respect to (x^, yr) = (4.0, 10.0) 

Case 2(c): {xp,yp) e i?2. For all {xp,yp) G i?2, the corresponding graph of -0 is 
or will eventually be bounded below by the graph of 0. Again, tolerance is possible 
but not guaranteed. 

In summary of initial condition (c) , given that the entire first quadrant is a region 
of inhibition, there is the possibility of tolerance for all {xr,yr) and {xp,yp) except 
when Xr < Xp < max(>o0i(i) and "0 is bounded above by 0, as illustrated in the 



orange region Ri in Figure p^c). Figure 3.6 links to three animations for Example 



3.17 ^c) with ^0(0) chosen from the different regions shown in Figure 3.4 'c). As be- 



fore, each animation shows phase space and time courses in a side-by-side comparison. 

We now use time interval estimates to expand the region that guarantees toler- 
ance. Consider initial value (a). We choose (x/,y/) such that j// — maxt>o02(i)- We 
note that the extremal points of 0(t), {xmtUm) and (xf,yf), are on the x-nuUcline 
and y-nuUcline respectively so that yM — xm — 1 and yf — 2a: f. Given that initial 



value (a) results in an excitable trajectory, we can apply (3.8 1 with Cr = 12 and 



= Cf = |a;^/(l-|- 2a: J ) — Xf\. This then establishes a bound on xm, such that 
tolerance occurs for Xr < Xp < xm, in terms of the initial value and extremal points 
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Fig. 3.6. Animations for Example \3 . li[ c) , showing the presence or absence of tolerance with 
respect to cj>(0) = (4, 10) for differing choices 0/1/1(0). </>(0) is denoted by the large red dot and ip{0) 
is denoted by the smaller blue dot. Given (f)(0) = (4,10), V(0) = (4.2,2) £ Ri does not produce 
tolerance (Left), '0(0) = (5,25) S R2 produces tolerance (Middle), and i/'(0) = (6,5) £ R2 does not 
produce tolerance (Right). 



of the reference trajectory (/'(i). For example, rough bounds on Xf and Xm can be 
obtained from a visual inspection of 4'(t). From Fig 3.4 we can propose 2 < Xf < 3, 
leading to 1.55 < C/ < 2.53, and 4.5 < xm < 5, with xm = xm + [C f /Cr){xM — Xr) 
from (3.8 1 with — Cf. More stringent bounds can be obtained by performing nu- 



merical integration using interval arithmetic. Moreover, as yp increases, increases 
while Cf remains fixed, such that tolerance can be guaranteed for larger Xp, given 
larger yp. 

In fact, example |3.17| is simple enough that we can obtain more precise estimates 
and t^, as defined in Section 3.3 Let (similarly, t^) be the time of passage 
(V'l = Xp) to (j)i — Xf ("01 — Xf). (j){t) can be represented by two 



on 
from 



with xj, 



= 1, 2 on the two 
Xf, where xm = 



for t e [O,i0] by {u,Vi{u)), i 
Xri x^ = Xm and = 
Recall that m this example, the entire first quadrant is a region of 
inhibition. Our approach is to estimate the time intervals by setting Vi{u) to a constant 
in (3.5 1 and then integrating to obtain > A(?/r, a;^, xm) + ^{yf,XM,Xf), where 



segments. Denote the gr aph of 
segments, is given by (3.5 1 
maxt>o0i(i) 



A(w, a, b) 



du 



\l + w — b\ , a /„ , , N 

u \l + w — a\ 



Next, we compute for the trajectory il'{t) with initial condition {xp,yp) and 
ending at {xf,yf). Now, consider those {xp,yp) such that Xp > Xr an d y^, > yf. Since 



the y-nuUcline is the curve y = 2x'^, by uniqueness of solutions to (3.9 1, the latter 
condition ensures that V'2(i) > j// for all t such that ipi{t) > Xf. By the continuity 
of A{w,Xf,Xp) in w, = A{y.^,Xp,Xf) for some y^ > yf. Thus, for the tolerance 
condition < to hold, it is sufficient that 



A{yr,Xr,XM) + A{y f , X M , X f) > A{y^,Xp,Xf). 



(3.12) 



If Xp — Xm, then the observation that A(ijf,XM,Xf) > A{y^, xm , Xf) implies that 
(3.121 holds, and hence tolerance occurs, 

Xp, Xf) 



> XM, writing A(y^, 



as expected from Proposition 3.10 
A(^y.^,Xp,XM) 



For 



A{y^,XM,Xf) shows immediately 



that the upper bound for tolerance can be extended from xm to some Xp > xm- 

Assuming that both sides are positive, as in Figure 3.4 condition (3.121 can be 
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expressed as 



Xr{l -Xf+ Vf){xM - 1 



(1 - XM + yf){Xr -l-Vr) 



Vr) ^ (1 + - Xf)Xp 



(3.13) 



Condition (3.131 still depends on y^, which can be estimated under the assumption 



that > ip2{t^) (which holds, for example, if g < along ■!/)(/:) from t = to 



t = t^). Formally integrating the second equation of (3.9 1 gives 4'2{t%ji) — ypS 



/g* e"*^**~*''/^a;^di'. On the trajectory ^{t), xj < x < Xp, hence ^/'2(iv) > 2/pe"*'''/^ + 
J^i- g^-(ti,~t )/2^2^^/ _ _|_ (^y^ _ y^-^^~t^/2^ whei B wc h a ve used y/ = 2x'j. Now 
tji, = A{y^,Xp,Xf) < A{yf,Xp,Xf). Therefore, > ^2{t^) > yb, where 

yb = yf + (yp - yf)exp[-A{yf,Xp,Xf)/2], 



(3.14) 



and t/h is an affine function of yp. Note that the right hand side of (3.13 1 is a monotonic 
decreasing function of y^j,. Hence, (3.131 is guaranteed to hold if 



Xr{l - Xf + yf){xM - 1 - 2/r) ^ (1 + J/b - Xf)Xp 



(1 - Xm + yf){Xr - 1 - yr) 



1 + yfc 



(3.15) 



which is a condition on tolerance for the initial value {xp,yp) of ip{t) in terms of the 
initial value and extremal points of the reference trajectory (/)(i). Finally, we note 



that condition (3.151 is also applicable for initial condition (b) or (c). In those cases, 
set Xm = Xr- 

Remark 8. Ifyp is increased for fixed Xp, then y^ increases, such that the right 
hand side of (3.13) decreases. Thus, the larger yp is, the more likely it is that (3.13) 
is satisfied. 

Example 3.18. Let y = rx—y, r > and consider the following general equations 
as possibilities for x — f(x, y): 



X = fix,y) 
X = f{x,y) 
X = f{x,y) = 



1 + by 
X = f{x, y)^ ax- by" 
ax" 



cx 



1 + by'' 



cx, 



(3.16) 
(3.17) 
(3.18) 



where a,b, c > 0, n, m e Z+ and x,y > 0. 



Each of the above equations models inhibition of x by j/, with /y < in the first 
quadrant, implying that the entire first quadrant is a region of inhibition. Assuming 
parameters are chosen so that (0,0) is a stable fixed point, results will be completely 
analogous to those in Example |3.17[ 

More diverse possibilities arise when fy > on at least a subset of the first 
quadrant. For example, suppose that f{x,y) is the product of two inhibitory terms, 
such as 



f{x,y) = {ax + by){ 



cx 



l + dy 



h), 
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with 6 < and a,c,d > 0. Indeed, 

sgn{fy) — sgn{cx{b — adx) + bh{l + dy)^). 

If /i > 0, then fy < for all {x,y) E M."^^, as in the previous example. If, however, 
h < 0, then fy changes signs in M.'^^. 

Example 3.19. Consider the nonlinear system 



X = f{x,y) = (.5a:-2/)(f^-l) 
y = 9{x,y) = OAx-y, 



(3.19) 



with (0,0) as a stable node. The isoclines for this system are the family of curves 
given by the equations: 

O 11 1 

^ —X - - + -C+ — v/4a;2 + 20a; + 30a;C + 25 + 50C+25C2, (3.20) 

O 11 1 

= — a; - - + -C- — ^4x2 + 20a; + 30a;C + 25 + 50C + 25C2, (3.21) 

Li Li i U 



where C € R. In Figure |3.7| the isoclines are drawn for various values of C € 
[—2.5, 5.0], in increments of 0.25. For each C S M, the two curves defined by equations 



(3.201 and (3.21 ) together form a continuous curve. A thick black curve in the figure 
emphasizes the two parts, with equation ( 3.20[ ) forming the curves above and equation 



(3.211 forming those beneath. The equation of this curve, which looks linear in the 
first quadrant, is given by y = -1 + .22'iQ^Jx{x + 2). 

The portion of the first quadrant containing the top portions of the isoclines is not 
a region of inhibition, since for fixed x, f is an increasing function of y there. However, 
the portion of the first quadrant containing the bottom portions of the isoclines is a 
region of inhibition, since / is a decreasing function of y there. The curves given by 
the portion of the x-nuUcline (C = 0) in the first quadrant are marked (red) to help 
delineate where the speed of the isoclines (i.e. x) is positive or negative. 



Figure 3.8 shows a specific solution, (/)(i), that will be considered for this example. 
The following features appear in Figure [3^ 

• (j){t) is the curve shown in black for initial condition 0(0) = (x^, y, ) = (2, 0.5). 

• The orange curve, denoted as (f), is the curve of points obtained by integrating 
(p{t) backwards in time from t = to t « —0.85, at which time it intersects 
the X-axis at i; w 2.5. 

• Let i?2 be the region shown in green together with the boundaries made by 
(1) the line segment {{x,y)\x = 2,0 < y < 0.5}, (2) the orange curve, (j), and 
(3) the X-axis. 

• Define the region Ri to be the complement of R2 in F^^ , namely Ri'^^T^^ \^^2- 

Recall that every point (xp,yp) will lie on or to the right of the line x — Xr, by 
{A3). The regions Ri and R2 are formed so that for {xp, yp) G V will be bounded 
below by (j) and for {xp, yp) £ i?2 \ </>, "0 will be bounded above by (j). The graph of 4>, 
in orange, creates a natural boundary (by uniqueness of solutions) between different 
classes of solutions ip{t). 

Case 1: Let {xp,yp) G -Ri. Then, t/j will be bounded below by 0. Note that the 
graph of (j) never enters the region of inhibition. Thus, any (xp, yp) £ Ri does not 
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C <0 




C<0 



C = 



c >0 



Fig. 3.7. Isoclines for Example \3. 19\ drawn for various values of Ci £ (—2.5, 5), in increments 
of 0.25. The thick black line marks the boundary of the region of inhibition. 



produce tolerance with respect to {xr,yr) by Proposition 3.14 



Case 2: Let (xp, i/p) g -R2 \ 0- The resulting t/i will be bounded above by cj). Thus, 
from Proposition 3.16 since there are no regions of inhibition that contain both ili{t) 
and for all t>Q, tolerance may occur for (xp,yp). However, if {xp,yp) lies on 
the orange curve <j), then ip{t) and are subsets of the same larger solution curve 
of the vector field (3.19) and both (t)i{t) and V'i(i) ~^ monotonically a,s t ^ 00. 
By Proposition 



3.1 



therefore, {xp,yp) will not produce tolerance. In addition, by 
continuity, there exists an open ball, B, around each {xp,yp) S 0, such that {xb,yb) 
will not produce tolerance for all {xf,,yi,) € B. Thus, the set of points which might 
produce tolerance is a strict subset of region R2 . This set can be characterized more 
extensively by two different arguments. 

First, it is clear that tolerance occurs if {xp,yp) = (2.0,0), since /(2,0) < 
/(2,0.5) < (see the animation associated with the middle panel of Figure 
Thus, tolerance occurs for all {xp,yp) in a ball around (2,0), intersected with F^p"^ . 
The speed f{xp, 0) becomes monotonically less negative as Xp increases toward 2.5, 



and tolerance does not occur for {xp,yp) = (2.5,0) by Proposition 3.1 Thus, toler- 



ance occurs for {xp,Q) for all Xp G [2, Xp) for some Xp G (2,2.5). Similarly, /(2, j/p) 
becomes monotonically less negative as t/p increases from 0, where tolerance occurs, 
to 0.5, where it does not. Hence, tolerance occurs for (2,yp) for all yp £ [0,yp) for 
some j/p G (0, 0.5). Therefore, there is a continuous curve connecting {xp, 0) to (2, j/p), 
call it Ct, such that tolerance occurs exactly when {xp,yp) is in the interior of the 
region bounded by {x — 2], {y = 0}, and Ct- 

Second, to definitively establish that tolerance occurs for some specific {xp,yp) G 
i?2, time interval estimates for specific trajectories must be made, as done in Example 
|3.17[ Figure |3.9| provides links to three animations for Example |3.19| using various 
choices of ip(0) from the different regions shown in Figure 3.8 



20 



J. DAY, J. RUBIN AND C.C. CHOW 




Example 3.20. Consider the nonlinear system: 

_ 1 Q I I 

(3.22) 



X = f{x,y) = a: (^^-1.9) 

y = 9{x,y) = x-y 



The isoclines for this system are the family of curves given by the equations 
.5 (19a: + IOC + V37x'^ - MOxC - 300C2) 



_ .5 (19a; + IOC - V37x^ - 340xC - 300C2) 

y = 9.T + IOC ' ^^-^^^ 
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Fig. 3.9. Animations for Example \3.19\ showing the presence or absence of tolerance with 
respect to 4>{0) = (2,0.5) for differing choices o/i/)(0). 0(0) is denoted by the large red dot and 
^l>{Q) is denoted by the smaller blue dot. Given 4i{0) = (2,0.5), ip{0) = (5,1) £ i?i does not produce 
tolerance (Left), V(0) = (2,0) S R2 produces tolerance (Middle), and tp{0) = (2.2,0.2) S R2 does 
not produce tolerance (Right). 



where C E R. In Figure [3.10| the isoclines are drawn in increments of 0.1 for values 
of C e [—1.2, 0] and in increments of 0.01 for C G [0, 1]. For C £ [0, 1], the two curves 



defined by equations (3.23) and (3.241 together form a continuous curve. The black 
line, y = 1, in the figure emphasizes the two parts, with equation (3.231 forming the 



curves above and equation (3.241 forming those beneath 



Basin Of Attraction 
for(0,0)in Q1 



Stable manifold of 
Saddle Point 




C=-1.2 



C<0 



C=0 

.C>0 

c=o 



C<0 



Fig. 3.10. Isoclines for Example \3.2(^ drawn for various values of C Si (—1.2, 1). 

A saddle exists at (0.72,0.72). The stable manifold of this saddle point forms a 
boundary for the basin of attraction of (0, 0), r(o.o). The blue shaded region in Figure 
3.10 shows the subset of r(o,o) in the first quadrant. A third fixed point (stable spiral, 
not labeled) in the first quadrant is located at (1.4,1.4), outside of Fi-n ni- The x- 



(0,0)- 

nuUclines (C = 0) are marked (red) to help delineate where the speeds associated 
with the isoclines (i.e. x) are positive or negative. 
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We define several disjoint subregions (see Figure 3.101 of the basin of attraction 
of (0, 0) in the first quadrant, as follows: 

• a - above (and including) the top component of the C = isocHne, 

• (3 - below the top component of the C = isocline and above (and including) 
the line y = 1, 

• 7 - below the line y = 1 and above (and including) the bottom component of 
the C = isocline, and 

• (5 - below the bottom component of the C = isocline. 



These subregions are relevant because C varies nonmonotonically in y for this 
example and are defined to assist with identifying regions of inhibition. If looked at 
separately, subregions a and /3 are both regions of inhibition and subregions 7 and 6 
are not regions of inhibition. However, additional complications may arise if (f> and ip 
are not in the same subregion on some time interval. 



C=-1.2 




(Xr^r)=(0.5, 0.5) 



Fig. 3.11. Top: For Example 3. 20\ possible {xp,yp) points fall in one of three regions: R\, the 
red area plus its boundaries, R2, the magenta region, and R3, the yellow region. 



Figure 3.11 shows one specific solution, (j){t) with 0(0) = (0.5,0.5), that will be 



considered for this example. The following features are also a part of Figure [3. 11| 

• The orange curve, denoted as (f>, is the curve of points obtained by integrating 
(j){t) backwards in time from i = to i w —1.75, at which time it intersects 
the X-axis at x w 1.0. 

• i?i is the region shown in red together with the boundaries made by (1) the 
line segment {{x,y)\x = 0.5,0 < y < 0.5}, (2) the orange curve (f>, and (3) 
the X-axis. 

To,o)\("Ui?i). 



i?2 is the region shown in magenta, defined as R2 =^ Ff 



• i?3 is the region shown in yellow to be R3 q) ^ ct- 

As usual, we consider points (xp, yp) that lie on or to the right of the line {x = x^}. 
The region Ri is formed so that for (xp, yp) di Ri \ (f), ip will be bounded above by (j). 
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For (xp,yp) G i?2 U R3, ip will be bounded below by (p. The graph of (p, in orange, 
creates a natural boundary (by uniqueness of solutions) for ip{t), as in the previous 
example. 



Case 1: Let {xp,yp) E Ri \ (/>. Then, ip will be bounded above by (p. Thus, 



from Proposition 3.16[ since there are no regions of inhibition that contain both ipit) 
and (f){t) for all t > 0, (xp,yp) might produce tolerance. This case is very similar to 
that considered in the previous example. Indeed, it is clear that tolerance occurs if 
ixp,yp) = (0.5,0), while tolerance does not occur if (xp,yp) lies on 0, by Proposition 
3.1 Again, there will be a continuous curve connecting {{x, y) : x = 0.5, < y < 0.5} 
to {{x,y) : 0.5 < a; < l,y = 0} such that tolerance occurs for all {x,p) in i?i below 
this curve and does not occur in Ri above this curve. Time interval estimates are nec- 
essary to prove that tolerance occurs or does not occur for specific choices of (xp, yp) 
in Ri. 



Case 2: Let {xp,yp) G i?2. Then, cp C S and tp will be bounded below by (p. Note 
that 7 is not a region of inhibition and that /3, although a region of inhibition by 
itself, has / > 0, such that no tolerance can occur before ?/' enters 6. But S is not 
a region of inhibition, and hence from Proposition 
produce tolerance with respect to {xr,yr)- 



3.14 



any {xp,yp) G R2 does not 



Case 3: Let {xp,yp) G R3. Since C < in i?3, it is possible in this case that 
tolerance will occur before ip leaves a. Alternatively, suppose that this does not 
happen. After ^p leaves a, it enters /3,7, and finally (5 as it converges toward (0,0). 
In theory, tolerance could occur after -0 enters S. However, is bounded below by 
(p and 6 is not a region of inhibition. Hence, as in Case 2, Proposition 3.14 implies 
that tolerance will not occur. In summary, if 0(0) G S and (xp,yp) G R3, then either 
tolerance occurs before ip leaves a or it does not occur at all. 




C=-1.2 



Fig. 3.12. Nonmonotonic convergence to (0,0) in Example \3^ 
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Fig. 3.13. Animations for Example \3.20\ using an alternative choice for 4>{0) , namely in region 
a, as shown in Figure \3.12\ The presence or absence of tolerance with respect to (f>{0) = (0.5, 2) 
is shown for differing choices of ip{0) in Region Rz- (piO) is denoted by the large red dot and V'(O) 
is denoted by the smaller blue dot. Given (l>{0) = (0.5,2), ^/'(O) = (0.7,4) S R-j produces tolerance 
{Left) and ip{0) = (0.7,3) G Rs does not produces tolerance (Right). 



Using the same nonlinear system given by ( 3.22 ) , consider an alternative choice for 
{xr,yr), namely one in a. Such a choice demonstrates some additional complexities 
that can arise in this type of example. Now, passes through regions where / < 
0, then / > 0, and finally / < again as it converges to (0,0). For different ip 
trajectories, either bounded above or below by (f) (see Figure 3.12), there are different 
time intervals when tolerance cannot occur or might possibly occur, which can be 
inferred from the isoclines. 

In the particular example shown, for the ip that is bounded below by 0, toler- 
ance cannot be ruled out in any region. In particular, let xm denote the x-value 
where (j) intersects the x-nullcline branch that forms the boundary between 7 and 5. 
If V'i(*) < when (j) passes from a to /3, then tolerance is guaranteed to occur. 
On the other hand, for the i/j that is bounded above by (j), tolerance is only possible 
after tp enters S. Figure [3 . 1 3 1 provides links for two animations for Example |3. 20| using 
4>{0) = (0.5,2) in Region and two choices of "0(0) also in Region R^, similar to 
those shown in Figure |3.12[ 



4. Tolerance in Linear ODE systems. The previous sections have established 
that it is sometimes difficult to make precise general statements about tolerance. 
However, in the case of linear systems, we can fully characterize the occurrence of 



tolerance for equation (2.1 1. In this section, we derive a complete set of necessary and 
sufficient conditions for the existence of tolerance in 2D linear ODE systems. 
Consider the linear system 

X = Ax, (4.1) 

where A G M^^^, x E = [0, 00) x [0, 00). Throughout this section, we will assume 
as before that: 
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(Al) (0,0) is a stable fixed point of (4.11, the eigenvalues of which are real and 
negative. 

(A2) (j)(t) and ipit) ^re nonnegative for all t > and both {xr,yr) and {xp,yp) lie in 

the basin of attraction for (0,0) in the first quadrant, FJj p^. 
(A.3) Xp ^ X'p. 

Let Al and A2 be the real, negative eigenvalues of A. To arrive at necessary and 
sufiicient conditions for the existence of tolerance, there are two cases that must be 
considered. The first case is that A has distinct eigenvalues, Ai 7^ A2. The other case 
is that A has identical eigenvalues, Ai = A2 li^ A < 0. For each of these cases, there 
are subcases to consider as well. 

4.1. Case 1: Ai 7^ A2. For this case, where Ai and A2 are distinct, negative 
eigenvalues of A, assume without loss of generality that A2 < Ai < 0. Let v be an 
eigenvector corresponding to Ai, and let w be an eigenvector corresponding to A2. 
Since Ai and A2 are distinct, v and w are linearly independent. Thus, any initial 
condition can be uniquely written as a linear combination of v and w. In particular, 
(xr, Ur) — civ + C2W — {ciVi + C2W1, C1V2 + C2W2), with Ci , C2 £ M. Then, the solution 
4){t) to the initial value problem (IVP) x = Ax, ^(0) — {xr,yr) is 

(pit) = ciwe^i* + C2we^'* = (ciwie^^* + C2Wie^^*, ciV2e^''^ + C2W2e^'*). (4.2) 

Similarly, consider the initial condition {xp,yp), which can be uniquely written as 
(xp, Vp) = div + = {divi + d2Wi, diV2 + ^21^2), with di,d2 G M. The solution ip{t) 
to the IVP X = Ax, "0(0) = {xp, yp) is 

V'(t) = diwe^i* + d2we^^* = (diwie^^* + d2Wie^-'\ diV2e^^^ + d2t«2e^"*). (4.3) 

Since we know Xp > Xr by (A3), we have that 

diVi + ^2^1 > cxvi + C2'Wx. (4.4) 

We will consider three subcases for Case 1: (a) ui = and wi = 1 (b) vi — \ and 
Wi = Q and (c) vi = wi = \. 



4.1.1. Case la: vi — Q and wi = 1. For this case, (4.4 1 becomes 

d2 > C2. (4.5) 



Consider the difference between (j)i(t) and ipiit)- Using equations (4.2 1 and (4.3 1 as 
well as = and wi = 1, we have 



By (4.51, we have that (c2 — ^2) < 0. Thus, because e'*'^* > for all i > 0, we have 
that 01 (t) — ipi{t) < for all t > 0. Therefore, the following result has been shown. 

Proposition 4.1. Assume (Al), {A2), {A3) and that X2 < \i < 0. Given 
{{xr,yr), {xp,yp)) , if vi = and wi = 1 for eigenvectors v and w of Ai and X2, 
respectively, then x = Ax does not exhibit tolerance for {(xr,yr), {xp,yp))- 



4.1.2. Case lb: vi = 1 and wi = 0. For this second subcase of Case 1, (4.4| 
becomes 



di > Ci. 



(4.6) 
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Using equations (4.2 1 and ( |4.3| as well as wi = 1 and wi — 0, we have 

0i(i)-Vi(i)-(ci-di)e^i*. 



By (4.6 1, we have that (ci — di) < 0. Thus, we conclude (f>i{t) — ipi{t) ^ for all 
t > 0, and the following result has been shown. 

Proposition 4.2. Assume (Al), {A2), {A3) and that A2 < Ai < 0. Given 
{{xr,yr),{xp,yp)) , if vi = 1 and wi = Q for eigenvectors v and w of Xi and A2, 
respectively, then x = Ax does not exhibit tolerance for {{xr, yr), {xp, j/p)). 

4.1.3. Case Ic: vi — wi — 1. Unlike Cases la and lb, tolerance is a possibility 
in case Ic. Proposition |4.3| below states necessary and sufficient conditions on the 
coefficients of the solutions (p and tp in order for tolerance to be exhibited and also 
specifies the precise time value beyond which tolerance is exhibited, when it occurs. 

Proposition 4.3. Assume (Al), {A2), {A3), A2 < Ai < 0, and wi = wi = 1 for 
eigenvectors v and w of Xi and X2, respectively. Given {{xr,yr),{xp,yp)) , then there 
exists T > such that (4-.1) will exhibit tolerance for all t > T if and only if Ci > di 
and C2 < d2- Furthermore, 



T = 



H{d2 - C2)/{C1 - dl)] 

Al — A2 



(4.7) 



Proof. Necessary Conditions. Assume that ci < di. Since vi 



Wl 



rewrite (4.4) as 



dl + ^2 > Ci + C2. 



1, we may 



(4.8) 



Consider the difference between (/'i(t) and tpi{t). Using (4.2), (4.3), vi = wi — 1, and 



(4.8), we have 



(t>i{t) - Mt) = (ci - di)e^^* + (C2 - d2)e^^* 
< (ci -di)e^i* + (di -ci)e^^* 
= (ci-di)(e^i*-e^^*). 

Since A2 < Ai < and ci < c?i, it follows that (f)i{t) — V'i(^) < 0, which means 
that ipiit) ^ 't'lit) for all t > 0. Hence, tolerance cannot be exhibited for ci < di. 



Similarly, it can be shown that (4.1) cannot exhibit tolerance for C2 > d2. Thus, 



Ci > di and C2 < ^2 are both necessary conditions for tolerance. 



Sufficient Conditions. Assume that ci > di and C2 < ^2 both hold. Using (4.2 1 



(4.3), and vi = wi = 1, we have 

Mt) - Mt) = (ci - di)e^^* + (C2 - d2)e^^* 



e^^'(ci - di). 



ci - di 

By assumption, (ci — di) > and (c2 — ^2) < 0, and thus 

e^^*(ci - di) > and ~ t! < 0- 
(ci - di) 
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Therefore, 



V (ci-rfi)/ 



(ci - rfl) 
ln[(rf2 - C2)/(C1 - di)] 



(Ai - A2) 

□ 

4.2. Case 2: Ai = A2 = A < 0. In this case, A has cither a one- or two- 
dimensional eigenspace. Thus, two subcases need to be considered. 

4.2.1. Case 2a: Ai = A2 = A < and A has a two-dimensional eigenspace. 

For this case, A is an eigenvalue of A with multiplicity two for which two linearly in- 
dependent eigenvectors can be found. Let v and w be linear independent eigenvectors 
of A. Then, any initial condition can be uniquely written as a linear combination 
of V and w. For the initial condition {xr,yr), we may write {xr,yr) = civ + C2W = 
{ciVi + C2Wi,CiV2 + C2W2), with ci,C2 G K. Thus, the solution, to the IVP 

X — Ax, (/)(0) = [Xr, Ur) is 

<j){t) = ciwe^* -I- C2we^* = (ciwi + C2Wi,ciV2 + C2W2)e^* = (2:^6^*, y^e^*). (4.9) 

Similarly, consider the initial condition (xp,yp), which may also uniquely be written 

as {xp,yp) = div + d2W = {diVi + d2Wi,diV2 + ^2^2), with di,d2 G K. The solution 
V'(t) to the IVP X — Ax, -0(0) — [xp, yp) is 

i}{t) = dive^* + d2we^^ = (diWi ^2^1, ^1^2 + ^2^2)6^* = (xpe^*, ?/pe^*). (4.10) 



Consider the difference between (j)i{t) and i/'i(^)- Using (4.9 1, (4.101, (A3) and 
the fact that e'^* > for alH > we have that : 

0l(t) - -01 (t) = Xre^* - XpG^* = {xr - Xp)e^^ < 0. 

Thus, tpi{t) > 4>i{t) for all < > 0, and the following has been shown: 

Proposition 4.4. Assume (AI), {A2), {A3) and that Ai = A2 = A < 0. Given 
{{xr, yr), {xp, yp)) , if X has two linearly independent eigenvectors, then x — Ax cannot 
exhibit tolerance for {xr,yr)- 

4.2.2. Case 2b: Ai = A2 = A < and A has a one-dimensional eigenspace. 

In this case, let v be an eigenvector of A. One solution to (4.11 is x^^^t) — ve^^ . A 
second solution to (4.1) is x'^'^^t) = vte^* + ue"^*, where u is a generalized eigenvector 
satisfying (A — XI)v — v. The initial condition (xr,yr) can be uniquely written as a 
linear combination of v and v, 

{Xr,yr) = ClV + C2V = {ciVi + C2Vi,CiV2 + C2W2),with Ci , C2 € M. 

The solution (j){t) to the IVP x — Ax, 0(0) — {xr,yr) is 

0(t) = ciue^* -I- C2{vte^^ + we^*) 

= (ciuie^* -I- C2(wite^* -t- wie^*), ciW2e^* + C2(w2ie^* + W2e^*))- (4.11) 
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Similiary, the initial condition, (xp, i/p), can be uniquely written as a linear combina- 
tion of V and v, 

{xp, Up) = div + d2V = {divi + d2Vi,diV2 + (i2W2),witli di,d2 e M, 

and the solution V'(^) to the IVP x — Ax, V'(O) = {xp,yp) is 

il^it) = dive^^ + d2{vte^^ + we^*) 

= {divie^' + d2{vM' + vie^'),diV2e^' + d2{v2te^' + ^26^*)). (4.12) 

The following proposition, given without the details of its proof, states the result 
for this case. 

Proposition 4.5. Let {{xr,yr),{^piyp))- Assume {Al), {A2), {A3), and that 
Ai = A2 = A < 0. Suppose that A has a one-dimensional eigenspace. Let v be an 
eigenvector of A and let v be a corresponding generalized eigenvector. 



(i) Lfvi = 1 and vi = 0, then there exists T > such that (4--1) will exhibit tolerance 
for {{xr, Ur), {xp, Up)) for all t > T if and only if ci < di and 02 > ^2 both 
hold. Furthermore, 

T='^, (4.13) 

C2 - d2 

and the difference between (j)i{t) ci^d '4'i{t) at t > T will be less than or equal 
to (c2 — d2)te^*. Therefore, m^ {(c2 — ^2)^6'^*} = ^^x?"' ''J^hich occurs at 

t = is the greatest degree of tolerance that is possible. 



(ii) Lf vi ^ 0, then (4--1) will not exhibit tolerance for {(xr,yr),ixp,yp)) . 



4.3. Eigenvector Configurations and Regions of Tolerance. Of the cases 
discussed above, only Cases Ic and 2b yield the possibility of tolerance. The results 
stated above give analytical conditions for the existence of tolerance in terms of coef- 



ficients of general solutions to (2.1 1. We find that these results are more useful when 
they are recast geometrically. To achieve this reformulation, we consider eigenvec- 
tor configurations (EVC) that accommodate solutions that satisfy the nonnegativity 



requirement (A2). Each such configuration is displayed in Figure 4.1 For each con- 
figuration, we subdivide the positive quadrant into regions and then, for {xr,yr) in 
each region, determine precisely which locations for {xp, yp) will lead to tolerance and 
which will not. The results for all the eigenvector configurations shown in |4.1| are 
summarized in Table [4T] and are illustrated in the figures referenced in the table. 



4.3.1. Eigenvector Configuration (a). For eigenvector configuration (a), seen 
in the top left panel of Figure |4.1[ there are three regions in which to consider initial 
conditions: 



• REGION la: {xr,yr) on the x-axis 

• REGION 2a: {xr,yr) in the first quadrant below the weak eigenvector v 
and above the x-axis 

• REGION 3a: {xr,yr) in the first quadrant above the eigenvector v 
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Fig. 4.1. Regions of interest in the first quadrant for four relevant eigenvector configurations. 
Note that we label the weak eigenvector v with one arrow and the strong eigenvector w with two 
arrows. 

Now, we explain how to identify the regions of tolerance given an initial condition 
{xj-tUt), using Regions la and 2a as examples. 

REGION la: First, we look at the case when the initial condition is on the 
a:-axis. In the top left panel of Figure |4.2| an arbitrary point on the a;-axis is shown 
in the context of eigenvector configuration (a), with lines drawn (portions dashed), 
showing the addition of scalar multiples of the two eigenvectors to attain the point 
{xr,yr)- We refer to these lines as the Ci-line and C2-line. In this case, they divide 
the first quadrant into three different subregions, as shown in the top right panel of 
Figure [42] 

Recall that the P trajectory's initial condition (xp, yp) was expressed as (xp, j/p) = 
div + d2W. For all {xp, yp) in a given subregion, there is a corresponding relationship 
between di , c?2 and ci , C2 . Using this relationship, we determine if there exists a region 
where the criteria ci > di and C2 < (^2 of Proposition |4.3| and the initial condition 
criterion {xp > Xr) are all satisfied. For any {xp,yp) in such a region, tolerance will 
occur, while for {xp,yp) not in such a region, tolerance will not occur . 

In fact, for eigenvector configuration (a), if {xr,yr) is on the x-axis, then there 
are no subregions in the first quadrant where both di < ci and ^2 > C2 hold. In 
particular, in Iia, di < ci and ^2 < C2; in Ilia, di > c\ and di < C2; and in I Ilia, 
di > Ci and di > Ci. Thus, there exist no {Xp,yp) that produce tolerance. 



REGION 2a: Let {xr,yr) be in the first quadrant below the weak eigenvector 
V (but not on the x-axis) in eigenvector configuration (a). The middle left panel 
of Figure 4.2 shows an arbitrary point in this region, with lines drawn (portions 
dashed), showing the addition of the two eigenvectors to attain the point [x^, yr)- The 
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Eigenvector 
Configuration: 



If {Xr,yr) is 

in Region: 



Then, tolerance 
is produced by 
{xp,yp) in: 



Figure 
Reference: 



(a) Figure 4.1 i 



(b) Figure 



4.13 



(c) Figure 



4.1 



la 



2a 



3a 



lb 



2b 



3b 



Ic 



2c 



None 

Region W2a 
Region IVsa 

Region I if, 
Region hb 
Region ll^i, 

Region IVic 
Region W2c 



Figure |4.2| 
(top) 

Figure |4.2| 
(middle) 

Figure |4.2| 
(bottom) 



Figure |4.3| 
(top) 

Figure |4.3| 
(middle) 

Figure |4.3| 
(bottom) 



Figure |4.4| 
(top) 

Figure |4.4| 
(bottom) 



(d) Figure 4.1 i 



Id 



Region Wid 



Figure 4.5 



Table 4.1 

Summary of tolerance results for eigenvector configurations shown in Figure \4~l\ 



middle right panel of Figure |4 . 2 1 shows the four subregions formed in the first quadrant 
by the ci-line and C2-line. Note that region TV 2a only includes points satisfying 
X > Xr- In general, we follow the convention of truncating those subregions that 
satisfy Proposition 4.3 to ensure that (A3) is satisfied. 



4.3 



fail and 
> Xr and 



In this case, if {xp,yp) ^ IIV2a, then the conditions of Proposition 
tolerance will not occur. In contrast, for (xp,yp) G W2a, we have that ^ _ 
that di < ci and d2 > C2, such that all of the conditions of Proposition 4.3 hold. 
Hence, for eigenvector configuration (a), if (xr,yr) is in the first quadrant below the 
weak eigenvector v (but not on the x-axis) , then tolerance will be exhibited precisely 
for all {xp,yp) € IV2Q. 
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REGION 3a: Similarly to the case of Region 2a, the ci-line and C2-line partition 
the first quadrant into four subregions, as shown in Figure |4.2[ The conditions for 
tolerance only hold in subregion IVsq , which has been truncated to include only points 
satisfying x> Xr- 




Fig. 4.2. Left Side: Eigenvector configuration (a) with an arbitrary initial condition {xr,yr) 
labeled in Region la-3a. Right Side: The first quadrant partitioned into several different subregions 
by the ci-and C2-lines associated with the point (xr,yr) = civ + C2W lying in one of the initial regions 
la-Sa. 
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Fig. 4.3. Left Side: Eigenvector configuration (b) with an arbitrary initial condition {xr,yr) 
labeled in Regions lb-3b. Right Side: The first quadrant partitioned into several different regions by 
the ci-and C2-lines associated with the point (a:,-, yr) = civ + C2W lying in one of the initial regions 
lb-3b. 

4.3.2. Eigenvector Configuration {h). For eigenvector configuration (6), seen 
in the top right panel of Figure |4.1| there are three regions in which to consider initial 
conditions: 

• REGION lb: (a;^, j/c) on the a;-axis 

• REGION 2b: {xc, j/c) in the first quadrant below the weak eigenvector v 
and above the a;- axis 

• REGION 3b: {xc, J/c) in the first quadrant above the weak eigenvector v 
and below the strong eigenvector w. 

The results for each region are summarized in Table |4.1| and shown in Figure |4.3[ 

4.3.3. Eigenvector Configuration (c). For eigenvector configuration (c), seen 
in the bottom left panel of Figure |4.1| there are two regions in which to consider 
initial conditions: 
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Fig. 4.4. Left Side: Eigenvector configuration (c) with an arbitrary initial condition (xr,yr) 
labeled in Regions lc-2c. Right Side: The first quadrant partitioned into several different regions by 
the ci-and C2-lines associated with the point (xr, yr) = cif + C2W lying in one of the initial regions 
lc-2c. 



• REGION Ic: {xr,yr) in the first quadrant below the weak eigenvector v 
and above the strong eigenvector w 

• REGION 2c: {xr,yr) in the first quadrant above both eigenvectors 

Table |4.1| along with Figure |4.4| summarize the conclusions about tolerance for the 
regions in eigenvector configuration (c). 




Fig. 4.5. Left Panel: Eigenvector configuration (d) with an arbitrary initial condition (xr,yr) 
labeled in Region Id. Right Panel: The first quadrant of eigenvector configuration {d) partitioned 
into four subregions by the ci- and C2-lines associated with the point {xr,yr) = civ + C2V lying in 
Region Id. 
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4.3.4. Eigenvector Configuration (d). To finish our analysis, we examine 



eigenvector configuration (d), seen in the bottom right panel of Figure 4.1 There 
is only one region in which to consider initial conditions to explore the existence of 
tolerance. 

• REGION Id: {xr,yr) in the first quadrant above v 
The conclusion regarding tolerance for this case (Case 2b) was given by Proposition 
|4.5[ which shows that it is necessary and sufficient that vi = 0, ci < di, and C2 > d2 



for tolerance to be exhibited in (4.11. In the left panel of Figure |4.5| an arbitrary 
point in Region Id is shown in the context of eigenvector configuration (d), with lines 
drawn (portions dashed), showing the addition of scalar multiples of the eigenvector 
V and the generalized eigenvector, v, to attain the point {x^,y^). Since Vi = Q was 
assumed, the blue line along the y-axis represents v. 

The conditions ci < di and C2 > ^2 ^^re satisfied precisely for those {xp,yp) € 
TVid, the region labeled in the right panel of Figure 4.5 Moreover, Xp > Xr in this 



region as well. Hence, tolerance will be produced by any {xp, yp) S IVid, when {x^, yr) 
is in Region Id under eigenvector configuration (d). 

5. Discussion and Conclusions. Our consideration of tolerance serves as an 
example of how dynamical systems questions can arise from biological phenomena. 
We initiated our analysis of tolerance under assumptions representative of typical 
experimental preconditioning protocols used in the study of the acute inflammatory 
response [5j El [SI [12l [16] . However, in this paper, we present a generahzed analysis, 
allowing relatively general choices of initial conditions for the reference and perturbed 
trajectories, since the ideas of inhibition and tolerance, as we have defined them, are 
themselves quite general. The goal of this analysis is to use information about the ini- 
tial conditions of the R and P trajectories and the vector field to determine a priori if 
the associated trajectories will or will not exhibit tolerance. In tolerance experiments, 
by applying the challenge dose to the preconditioning trajectory at different times, 
an experimentalist could generate a continuous curve of possible initial conditions for 
what we call the P trajectory, and our analysis aims to consider all such initial con- 
ditions, to fully characterize the possibility of tolerance within a given experimental 
set-up. 

In the context of two-dimensional nonlinear systems of ODE, it can be difficult 
to malce general statements specifying conditions under which tolerance will be guar- 
anteed to occur. However, our work provides several fundamental statements about 
configurations of the initial condition (xp,yp) for the P trajectory, relative to the R 
trajectory, that will or will not lead to tolerance. For example, in Section [3T| we have 
characterized the case when the R trajectory is n-excitable, showing that there always 
exists a subset of the basin of attraction where tolerance is guaranteed to occur for all 
{xp,yp) in the subset. Excitable trajectories are common in systems describing vari- 
ous biological constructs and the idea of tolerance may be important to the ensuing 
analysis of such systems. By using isoclines and the concept of inhibition, we also 



present a framework in Section 3.2 that can be used to derive specific conditions under 
which tolerance can be ruled out or guaranteed in particular examples. Techniques 
such as time interval estimates in Section |3.3| exploit these ideas to achieve a closer 
examination of transient behavior in the absence of an analytical solution. 

In the linear case, we have fully characterized the conditions under which tolerance 
will or will not occur. A graphical view of the phase plane immediately reveals points 
'xp,yp) that produce tolerance relative to a given (xr,yr)- For example. Figures 



4.2- 



4.5 show regions of {xp,yp) (marked in green and labeled) in which tolerance will be 
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exhibited. Interestingly, some of the tolerance regions shown have infinite area (see 



Figures 4.3 4.4 and 4.51. Considering points {xp,yp) in the first quadrant and to 



the right of the vertical line x = x^, we see that in most cases (for instance, see the 
panels in Figure 4.3), the farther Xp is from Xr, the higher the Up value needs to be in 
order for (xp,yp) to fall in the green tolerance region. (As shown using time interval 
estimates this is also true in nonlinear systems.) Correspondingly, for some (xp, j/p) 
in a tolerance region, tolerance might only occur in the asymptotic limit, which may 
not be of interest in applications, especially considering that the degree or magnitude 
of tolerance produced is negligible by then. In other examples (for instance, see the 
middle and bottom panels of Figure 4.2), the value needs to be sufficiently low for 
tolerance to occur, although there is a limit on how low it can be because of the 
non-negativity requirement on y. 

The issue of tolerance, as defined in this work, does not appear to have received 
previous analytical treatment. Research has been done on isochronicity, which con- 
siders whether multiple phenomena occur within the same interval of time |10l |S] . For 
instance, in |10j , Sabatini defines a critical point classified as a center to be isochronous 
if every nontrivial cycle within a neighborhood of the critical point has the same pe- 
riod. Although Sabatini noted that the definition of isochronicity does not require 
proximity to a critical point, his work and other previous research appears to have 
been restricted to locating isochronous sections of autonomous differential systems 
that are oscillatory in nature ^\ |SJ [7]. While tolerance is a natural extension of 
isochronicity, in that it can be cast in terms of a comparison of the relative passage 
times of trajectories between sections, previous work has not, to our knowledge, made 
such comparisons between trajectories converging to a stable node, as we have done 
here. 

Another related area of study is the consideration of phase response curves (PRCs), 
as are commonly used in the analysis of neuronal systems. PRCs are calculated to 
determine how instantaneous perturbations shift the phase of a periodic oscillation. 
Although the assumption of intrinsic oscillatory behavior distinguishes the use of 
PRCs from the tolerance phenomenon that we consider, a relation between the two 
emerges if one thinks of an instantaneous perturbation as a preconditioning event 
and considers how the subsequent dynamics, during a specific window of time, com- 
pares to the unperturbed oscillation. Depending on where the perturbation occurs 
in the oscillation cycle, the occurrence of a stereotyped event, such as a peak, can 
be advanced or delayed relative to the unperturbed case, and the former could be 
considered as a form of tolerance, in that it would represent a speeding up of the 
event of interest. Figure [5?T] illustrates an example of such a phase advance, using the 
Morris-Lecar model. In theory, isoclines could be used to predict whether perturba- 
tions in a given system speed up or advance an oscillation. Past work has pointed 
out that PRCs corresponding to infinitesimal perturbations are intimately related to 
isochrons, or curves of constant asymptotic phase [T5], but these are different than 
isoclines. Indeed, analysis developed previously for PRCs (see e.g. [5] for a review) 
sheds little light on tolerance under the assumptions that we consider, since there 
is no intrinsic oscillation involved here. Note that the absence of an oscillation is 
quite characteristic of the types of models that motivated this work (e.g. [5]), since 
perturbations typically lead to a non-oscillatory decay to a healthy critical point or 
approach to one or more unhealthy, perhaps lethal, critical points. 

The work presented here looks exclusively at two dimensional ODE systems. Some 
of the results and techniques considered do not naturally extend to higher dimensions. 
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Fig. 5.1. Tolerance in the voltage equation of the Morris-Lecar model seen during a specific 
comparison time window. 

unfortunately. In [5] it was shown that the presence and magnitude of tolerance in 
a four dimensional ODE model of the acute inflammatory response depended not 
only on inhibition but also on the relative levels of the variable being inhibited when 
various doses of endotoxin were administered, through various feedback effects in the 
system. In the 2D linear case, the relationship between the level of the inhibitory 
variable and the relative level of the inhibited variable is most clearly seen. Refining 
the results for the 2D nonlinear case and extending the results for both linear and 
nonlinear systems to dimensions greater than two remains to be done. The present 
work, however, yields new and potentially useful insight into the behavior of transients 
away from the critical points to which they eventually converge, in the context of some 
types ODE systems that commonly arise in models of biological systems. 
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